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Introduction. 

1. This paper treats of the propagation of vibrations over the surface of a " semi- 
infinite " isotropic elastic solid, i.e., a solid bounded only by a plane. For purposes 
of description this plane may be conceived as horizontal, and the solid as lying below 
it, although gravity is not specially taken into account. # 

The vibrations are supposed due to an arbitrary application of force at a point. In 
the problem most fully discussed this force consists of an impulse applied vertically to 
the surface ; but some other cases, including that of an internal source of disturbance, 
are also (more briefly) considered. Owing to the complexity of the problem, it has 
been thought best to concentrate attention on the vibrations as they manifest 
themselves at the free surface. The modifications which the latter introduces into 
the character of the waves propagated into the interior of the solid are accordingly 
not examined minutely. 

The investigation may perhaps claim some interest on theoretical grounds, and 
also in relation to the phenomena of earthquakes. Writers on seismology have 
naturally endeavoured from time to time to interpret the phenomena, at all events in 
their broader features, by the light of elastic theory. Most of these attempts have 
been based on the laws of wave-propagation in an unlimited medium, as developed 
by Green and Stokes ; but Lord Rayleigh's discovery f of a special type of surface- 
waves has made it evident that the influence of the free surface in modifying the 
character of the vibrations is more definite and more serious than had been suspected. 
The present memoir seeks to take a further step in the adaptation of theory 
to actual conditions, by investigating cases of forced waves, and by abandoning 
(ultimately) the restriction to simple-harmonic vibrations. Although the circum- 
stances of actual earthquakes must differ greatly from the highly idealized state of 

* Professor Eromwich has shown ('Proe. Lond. Math. Soc.,' vol 30, -p. 98 (1898)) that in such 
problems as are here considered the effect of gravity is, from a practical point of view, unimportant. 
t 'Proc. Lond. Math. Soc./ vol. 17, p. 4 (1885) ; ' Scientific Papers/ vol. 2, p. 441. 
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things which we are obliged to assume as a basis of calculation, it is hoped that the 
solution of the problems here considered may not be altogether irrelevant. 

It is found that the surface disturbance produced by a single impulse of short 
duration may be analysed roughly into two parts, which we may distinguish as the 
" minor tremor" and the "main shock," respectively. The minor tremor sets in 
at any place, with some abruptness, after an interval equal to the time which 
a wave of longitudinal displacement would take to traverse the distance from the 
source c Except for certain marked features at the inception, and again (to a lesser 
extent) at an epoch corresponding to that of direct arrival of transversal waves, it 
may be described, in general terms, as consisting of a long undulation leading up to 
the main shock, and dying out gradually after this has passed. Its time-scale is 
more and more protracted, and its amplitude is more and more diminished, the 
greater the distance from the source. The main shock, on the other hand, is pro- 
pagated as a solitary wave (with one maximum and one minimum, in both the 
horizontal and vertical displacements) ; its time-scale is constant ; and its amplitude 
diminishes only in accordance with the usual law of annular divergence, so that its 
total energy is maintained undiminished. Its velocity is that of free Ilayleigh waves, 
and is accordingly somewhat less than that of waves of transversal displacement in 
an unlimited medium. # 

The paper includes a number of subsidiary results. It is convenient to begin by 
attacking the problems in their two-dimensional form, calculating (for instance) the 
effect of a pressure applied uniformly along a line of the surface. The interpretation 
of the results is then, comparatively simple, and it is found that a good deal of 
the analysis can be utilized afterwards for the three-dimensional cases. Again, the 
investigation of the effects of a simple-harmonic source of disturbance is a natural 
preliminary to that of a source varying according to an arbitrary law. 

Incidentally, new solutions are given of the well-known problems where a periodic 
force acts transversely to a line, or at a point, in an unlimited solid. These serve, to 
some extent, as tests of the analytical method, which presents some features of 
intricacy. 

2. A few preliminary formulae and conventions as to notation may be put in 
evidence at the outset, for convenience of reference. 

The usual notation of Bessel's Functions " of the first kind " is naturally adhered to ; 
thus we write : 

Jo (0 = ~-~ cos (£ cos w ) ^ w ........ (i). 

TT J() 

* Compare the concluding passage of Lord Rayleigh's paper : 

" It is not improbable that the surface-waves here investigated play an important part in earthquakes, 
and in the collision of elastic solids. Diverging in two dimensions only, they must acquire at a great 
distance from the source a continually increasing preponderance." 

The calculations indicate that the preponderance is much greater than would be inferred from a mere 
comparison of the ordinary laws of two-dimensional and three-dimensional divergence. 
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By a known theorem we have also 

J (£) = sin (£ cosh u) du (2), 

IT Jo 

provided £ be real and positive. For our present purpose it is convenient to follow 
H. Weber # in adopting as the standard function " of the second kind" 

2 f 00 
K (£) = — cos (£ cosh u) du (3). 

IT J() 

It is further necessary to have a special symbol for that combination of the two 
functions (2) and (3) which is appropriate to the representation of a diverging wave- 
system ; we write, after Lord R,AYLEiGH,f 

D (£) = 2 ( e~^ coshu du . ........ (4), 

rr Jo 

so that 

D o (0 = K o (£)-a o (O (5)- 

We shall also write, in accordance with the usual conventions, 

J 1 (O=-J / (a K l (C) = -K' (O, Di(Q=-VoU) ■ • (6). 
For large values of £ we have the asymptotic expansion 

DbW -Vs ■•■"*" I 1 -UMT + 2175*? — •}• ' ' <7) ' 



In the two-dimensional problems of this paper we shall have to deal with a 
number of solutions of the equation 

3x 3 3v 3 

constructed from the type 

where <f is real, and 

(f*-A*), or = V(A a -f 2 ) (10), 



a = v 

* ' Part. Diff.-Gleichungen d. math. Physik,' Brunswick, 1899-1901, vol. 1, p. 175. Heine (' Kugel- 
functionen,' Berlin, 1878-1881, vol. 1, p. 185) omits the factor 2/ir. In terms of the more usual 
notation, 

K = 2 {- Y + (log2-y)Jo}, 

where y is Euler's constant. The function |ttK has been tabulated (see J. H. Michell, ' Phil. Mag.,' 
Jan., 1898). 

t 'Phil. Mag./ vol. 43, p. 259 (1897); 'Scientific Papers,' vol 4, p. 283. I have introduced the 
factor 2/?r, and reversed the sign, 

B 2 



4 PROFESSOR HORACE LAMB ON THE PROPAGATION OF 

according as £ 3 ig h*> the radicals being taken positively. In particular, we shall 
meet with the solution 

, _ 1 r e" ay e^ x d^ _ 2 f 30 e~ fly co s gag c?£ , . # 

and it is important to recognize that this is identical with D (Ar) ? where 
r= ^/(x % + ^/ 3 ). To see this, we remark that <f> 9 as given by (11), is an even 
function of x 9 and that for x = it assumes the form 



2 



iOO 



^™ a2/ <*f _ 2 n e-^ (w+ ^^ 



by the method of contour-integration.* This is obviously equal to D (%), Again, 
the mean value of any function <f> which satisfies (8), taken round the circumference 
of a circle of radius r which does not enclose any singularities, is known to be equal 
to J (kr) t <f) , where <£ is the value at the centre, t We can therefore adapt an 
argument of Thomson and TaitJ to show that a solution of (8) which has no 
singularities in the region y > 0, and is symmetrical with respect to the axis of y, is 
determined by its values at points of this axis. We have, accordingly, 



D (hr) = X f £l^jk ........ (13). 

Again, in some three-dimensional problems where there is symmetry about the 
axis of z, we have to do with solutions of 

J _ ,«*Ls» I, n-la- „__jtL. «„L, /} ^ft% *~*"»=» ( ) / ] A § 

3 j— ^ J~ ^ —J— Ifj <UJ •__ \/ „ . , . , , , IX'Xfj 

based on the type 

$ = Ae""" ft2 J ($er) ......... (15), 

where ^ = v /(^ s + y s ), and a is defined as in (10). In particular, we have the 
solution 

fh _ r°°^ as Jo.(^)i ^| / . fi x 

J Q Ot 

which (again) reduces to a known function. At points on the axis of symmetry 
(ot = ()) it takes the form 

Jo GL lih Z 

* If we equate severally the real and imaginary parts in the second and third members of (12), we 
reproduce known results. 

f H. Weber, ' Math. Ann., 5 vol. 1 (1868), 
J c Natural Philosophy/ §498, 
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Since the mean value of a function <f) which satisfies (14), taken over the surface of 
a sphere of radius r not enclosing any singularities, is equal to 

sin hr , 
fir 

where <f> is the value at the centre, # the argument already borrowed from Thomson 
and Tait enables us to assert that 



e 



— ihr 



r 



.00 







e-"J (^j)ed$ 



a 



( | o\ 4- 

a • • • • • • • \ -*- ^ ' } I 



where 



V = v/(^ + s») ^ ,/(^ + ^ + Z *). 



Finally, we shall require Fourier's Theorem in the form 



1 



.30 



,"»00 



/(«•)= ~\ rff 



/(\) «*<*-*> tf\ 



9 « 9 



X «/ GO 



« * 



• (19)4 



and the analogous formula 



.00 



'GO 



/(*,) = J (Hf^ f(\)J Q (fr)\dk (20).§ 

'Jo J () 

As particular cases, if in (19) we have /(a?) = 1 for x l < a 2 , and = for x 2 > a 3 , then 



.-•CO 



/(*) 



7T 



Sill L0& ,;c,. 7 > 

-co g 



2 

7T J 



/•CO 







sin £a 



cos £x d£ . . . . (21) ; 



and, if in (20) /(sr) = 1 for m < a, and = for sy > a, then 



/■•CO 



/(«r) = a J (far) J x (£a) c££ . 

Jo 



• • . • . « \ZiZj). 



These are of course well-known results. 



* H. Weber, < Crelle,' vol. 69 (1868). 

t If in (18) we put z — 0, and then equate separately the real and imaginary parts, we deduce 

I Jo (£ cosh u) cosh u du = ■-'-— , 



00 







sin i 

J (C sin w) sin w du = — ^ . 



These are known results. C/. Rayleigh, ' Scientific Papers/ vol. 3, pp. 46, 98 (1888) ; Hobson, ' Proc. 
Lond. Math. Soc.,' vol. 25, p. 71 (1893) ; and Sonine, < Math. Ann.,' vol. 16). 

J H. Weber, 'Part. Diff.-Gl. etc.,' vol. 2, p. 190. Since A, occurs here and in (20) only as an inter- 
mediate variable, no confusion is likely to be caused by its subsequent use to denote an elastic constant. 

§ H. Weber, 'Part. Diff.-Gl. etc.,' vol. 1, p. 193. 

|| It may be noticed that if in (20) we put f(vr) = e-^^/w, we reproduce formulae given in the foot-note t 
above. 



6 



PROFESSOR HORACE LAMB ON THE PROPAGATION OF 



PAET I. 



m 



Two- Dimensional Problems. 



3. The equations of motion of an isotropic elastic solid in two dimensions (x y y) are 



f)%1l 3A cfti) 3A 

p-^2 = ( X + /*) *„ + f*V%, p - d -- 3 - = (X + ft) ^ + /^V 2 v . . (23), 



df 



dx 



dP 



where u, v are the component displacements, p is the density, X, /a are the elastic 
constants of Lame, and 

(24). 



du . dv 
dx dy 



9 



« 6 • 



These equations are satisfied by 



u 



provided 



3<£ , d\jj 
dx dy 



v 



d(b dxb 
dy 



* 



& 8 9 9 



X 



dt 2 



p at" 



/x 



V^XJJ . 



# a s a 



(25),* 



(26). 



In the case of simple-harmonic motion, the time-factor being e^, the latter 
equations take the forms 



(V 2 + A 2 ) (£ = 0, (V 2 + k") xp = 



s e s 3 6 



where 



/^ 



P 2 /> 



A -p 2/x 



p 2 a 2 . 



P^P 



pW . 



. (2/), 



(28), 



the symbols a, 6 denoting (as generally in this paper) the wsbve-slownesses^ i.e., the 

reciprocals of the wave-velocities, corresponding to the irrotational and equivoluminal 
types of disturbance respectively. 

The formula (25) now give, for the component stresses, 



Pxx 

P X JI 

(I 

Pl/l 

fJL 



X A + 2^ 

/X OX 

dv M du _ 
dx dy 

^A + 2? V 



¥■<!> 



9 &<f>_ _L 9 ^ ... ^ 

r\ T f\ f\ 

O'lf ox oy 



2 -^ 

dxdy 



Jvfy 



» 3 2 'l- 



3x 



(29). 



#ty - 2 



2 



3.x 



C 



2 . !¥ 

dxdy 



* Green, 'Camb. Trans./ vol. 6 (1838); * Math. Papers/ p. 261. 

t The introduction of special symbols for wave-slownesses rather than for wave-velocities is prompted by 
analytical considerations, The term " wave-slowness " is accredited in Optics by Sir W. B. Hamilton. 
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In the applications which we have in view, the vibrations of the solid are supposed 
due to prescribed forces acting at or near the plane y = 0. We therefore assume as 
a typical solution of (27), applicable to the region y > 0, 

<£ = Ae-te**, xjj = Be-toe** ....... (30), 

where £ is real, and a, ft are the positive real, or positive imaginary,^ quantities 

determined by 

a a = p _ h\ £ s = f* -**...... . (31). 



For the region y < 0, the corresponding assumption would be 

<\> = k!e* v S\ y\s = B'ePe*** 



. (32). 



The time-factor is here (and often in the sequel) temporarily omitted. 

The expressions (30), when substituted in (25) and (29), give for the displacements 

and stresses at the plane y = 

u () = (if A ~~ ftB) e«*, v = ( - «A - #B) e** .... (33), 
and 

Mo = /* {(21 s - A 8 ) A + 2i£8B} e* 



(34). 



The forms corresponding to (32) would be obtained by affixing accents to A and B, 
and reversing the signs of a. ft. 

4. In order to illustrate, and at the same time test, our method, it is convenient to 
begin with the solution of a known problem, viz., where a periodic force acts 
transversally on a line of matter, in an unlimited elastic solid, f 

Let us imagine, in the first instance, that an extraneous force of amount Ye*** per 
unit area acts parallel to y on a thin stratum coincident with the plane y = 0. The 
normal stress will then be discontinuous at this plane, viz., 



, Pyy 



y=-+o 



L Pyyly ~ 







Ye^ x 



(35), 



whilst the tangential stress is continuous. These conditions give, by (34). 

Y 



(2P - F) (A - A') + 2i£8 (B + B') = 



jU, 



2ifr (A + A') + (2^ 3 - **) (B - B') = 



• (3 6 )- 



Again, the continuity of u and v requires 

ig (A - A') - /3 (B + B') = 
«(A + A') -f^(B- B') = 



I. 



(37). 



* This convention should be carefully attended to ; it runs throughout the paper. 
t Eayleigh, ( Theory of Sound, 5 2nd ed., § 376. 
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Hence 

A — . A' — ^- "R — - TV — - -^ __ /^^ 

xx ___» — — _£j_ __ -- - j jl> ~—«- .O —"*< '1 • rjTO"' s * * " ° \ /" 

Zirv Lb £j ZifC Us 

We have, then, for y > 0, 



*£ 



To pass to the case of an extraneous force Q concentrated on the line # = 0, ?/ = 0, 
we make use of (19). Assuming that the /(A.) of this formula vanishes for all but 
infinitesimal values of X, for which it becomes infinite in such a way that 



>QO 



/ I A. J Cv/v «—»*>•< \cvj i 



•00 



we write, in (39), Y = Q dij/2ir, and integrate with respect to £ from — go to + qo , # 
We thus obtain, for y>0, 

A — 0. . [ p -ni f A* f ]£ ,/, — __?Q f &J^d% (4Q\ 

47T/C jLL J — oo 47T/U /X - -co p 

or, on reference to (13), 

^-1^A,(H ^=^ ^Do^O- • • • • (4i), 

where r = ^/(x 2 -j- i/ 2 ). 

If we put x = r cos 9, y = r sin #, we find from (25), on inserting the time-factor, 
that for large values of r the radial and transverse displacements are 



3i/r Q / 2 

r30 dr 4/x V ^ 



a <£ _ ^ _ Q A / 2 ^-ir-i.0 cos # 



* * 



(42), 



respectively, t Use has here been made of (7). 

A simple expression can be obtained for the rate (W, say) at which the extraneous 

* The incleterminateness of the formula (19) in this case may he evaded by supposing, in the first 
instance, that the force Q, instead of being concentrated on the line b = 0, is uniformly distributed over 
the portion of the plane # = lying between x= ±a. It appears from (21) that we should then have 

If we finally make a = Q we obtain the results (40). 

t The second of these results is equivalent to that given by Rayleigh, he, tit, for the case of 
in compressibility (A, = go ). 
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9 



force does work in generating the cylindrical waves which travel outwards from the 
source of disturbance. The formula) (40) give, for the value of dv/dt at the origin, 



dv _ ipQeP* 
dt Ank^/i 



J. 00 



J —00 



It 



— a)dg 



» • 



» 8 



. (43). 



This expression is really infinite, but we are only concerned with the part of it in the 
same phase with the force,* which is finite. Taking this alone, we have 



d u _ pQeP* | 
dt \ttH 1 ^ \ 



k 



i y(P-e) + L^ h '-^ dl h^+ h ^ ■ < 44 >- 



Discarding imaginary parts, we find that the mean rate, per unit length of the axis 
of z, at which a force Q cos pi does work is 

-._>+ 3 tf p QZ (45). 



W 



^WJ 16 [i 16/x(A. + 2/x) 



5. We may proceed to the case of a " semi-infinite ?? elastic solid, bounded (say) 
by the plane y = 0, and lying on the positive side of this plane. We examine, in 
the first place, the effect of given periodic forces applied to the boundary. 

As a typical distribution of normal force, we take 



-PwJQ — V, [_Pyy„ 







Y,Mx 
C- * * • . • 



e » « 



(46), 



the factor e'# being as usual understood. The constants A, B in (30) are determined 
by means of (34), viz. : 



Hence 



2ifrA + (2f 2 - F) B = 0, 

(2£ s - /c 2 ) A + 2t£8B = X 



s • 



~/ 



• . « 



• II 



(47). 



A 



where, for shortness, 



2g 8 - ¥ _ Y 



~F(0'//-' ' 



e « 



F (£) = (2| 2 - £ 2 ) 2 - 4£ 2 «/3 



»#*•«•$ 



(48), 



(49). 



We shall find it convenient, presently, to write also 



/(£) = (2P-F)2 + 4£ a a£ (50). 



The awkwardness is evaded if (as in a previous instance) we distribute the force uniformly over a length 
2a of the axis of x. This will introduce a factor /SH_l?y un dei* the integral signs in the second 



member of (44). 
VOL. coin. — A. 







TO 



PROFESSOR HORACE LAMB ON THE PROPAGATION OF 



The surface- values of the displacements are now given by (33), viz, : 



u 



o 



_ if (2p _ P - 2a/3) ^ Y 



•\ 



F (f ) 



V, 







F (£) fi 



r 



• • • I c/ JL ?» 



_> 



The effect of a concentrated force Q acting parallel to y at points of the line 
# = 0, ?/ = is deduced, as before, by writing Y = — Qd^/2ir 9 and integrating 

from - oo to oo : thus 



u 







27T/X J— "~ ^ ~F~(£J" 



~\ 



v o — 



Q 



„00 



hPae 1 ** dP 



r 



« • e \ *^ **^ / * 



2tt/x J-.0, F (£) 



J 



In a similar manner, corresponding to the tangential surface forces : 



we should find 






• (53), 






2i£§ X 

F(ir /*' 



2£ s - 



7 ", 
w 



~L 



fl 



. (54). 



And, for the effect of a concentrated force P acting parallel to x at the origin. 



u ^ 



v — 



2?r/x 
iP 



J — 00 



e<x> 



ZjiTLki *> 



F/3e fo dg 
f(2f 3 - 



1 



J- 



•X) 



F — 2_a/3) f^d,( 



. « I iJ «J Is 



The comparison of a in (52) with w in (55) gives an example of the general 
principle of reciprocity. # 

We may also consider the case of an internal source of disturbance, resident 
(say) in the line x = 0, y =f, the boundary y = being now entirely free. The 
simplest type of source is one which would produce symmetrical radial motion (in 
two dimensions) in an unlimited solid, say 



(f) = D (hr), xfs = 



» » 



• « 



(56) : 



where r, = <y/{ s2 + (?/ ~~" jQ 2 }> denotes distance from the source. If we superpose on 
this an equal source in the line x = 0, y = —f 9 we obtain 

tf> = I) (hr) + T> Q (hr'), $ = ...... . (57), 



* Kayleigh, 'Theory of Sound/ vol 1, § 108. 
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where r f = ^'{x 9, + (y + jQ 3 }. It is evident, without calculation, that the condition 
of zero tangential stress at the plane y = is already satisfied ; the normal stress, 
however, does not vanish. It appears from (13) that in the neighbourhood of the 
plane y = the preceding value of <f> is equivalent to 



-00 



OL IT J — od Ot 



f ^e-^*^ (5 8 ). 



Zi 

7T J — oo CC 



Substituting in (29) we find that this makes 



9// f 00 9£2 _ 7,2 

>Jo = 0, Mo = 4M -— -^e-^6«"^ .... (59). 



7r J—oo a 



Comparing with (46), we see that the desired condition of zero stress on the 
boundary will be fulfilled, provided we superpose on (57) the solution obtained from 
(30) and (48) by putting 



rr a 



and afterwards integrating with respect to f from — go to oo . The surface- 
displacements corresponding to this auxiliary solution are obtained from (51), and if 
we incorporate the part of u due to (58), we find, after a slight reduction, 



u 



% 



4* 2 r P&r* e** d( 

^ . . . . . . (^0)« 

2/c 3 r (2g 3 - k?) e~ af e^ dg 

7T J-oo F (f) y 



. These calculations might be greatly extended. For example, it would be easy, 
with the help of Art. 4, to work out the case where a vertical or a horizontal periodic 
force acts on an internal line parallel to z. And, by means of the reciprocal theorem 
already adverted to, we could infer the horizontal or vertical displacement at an 
internal point due to a given localized surface force. 

6. It remains to interpret, as far as possible, the definite integrals which occur in 
the expressions we have obtained. 

It is to be remarked, in the first place, that the integrals, as they stand, are to a 
certain extent indeterminate, owing to the vanishing of the function F (£) for certain 
real values of f. It is otherwise evident a priori that on a particular solution of any 
of our problems we can superpose a system of free surface waves having the wave- 
length proper to the imposed period 2ir/p. The theory of such waves has been given 

c 2 
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by Lord Rayleigh,* and is moreover necessarily contained implicitly in our 
analysis. 

Thus, if we put Y = in (47), we find that the conditions of zero surface-stress are 
satisfied,, provided 

A : B = 2/c 3 — k 2 : 2&Vca-, = — 2ik(3 1 : 2k 3 — k 2 . . . (61), 

where k is a root of F (£) = 0, and et l9 /3 l5 denote the corresponding values of a, /?, 
Now, in the notation of (49) and (50), 

f (f)/(f) = (2^ - P)* _ i6 (f ~ ^) (e - v) e 

= / £ »{l-8|+(24-I 6 |)f,-l 6 (l.-|)|}. . (62). 

Equating this to zero, we have a cubic in £ s //c 3 , and since & 2 > A 3 , it is plain that 
there is a real root between 1 and co . It may also be shown without much difficulty 
that the remaining roots, when real, lie between and h?jk % . The former root makes 
a, /3 real and positive, and therefore cannot make /(f) = 0. The latter roots make 
a, j3 positive imaginaries, and therefore cannot make F (£ ) — 0. This latter 
equation has accordingly only two real roots f = ± *? where k > k. 
Thus, in the case of incompressibility (X = oo , h = 0) it is found that 

zc/& i= 1*04678 . . . , 

and that the remaining roots of (62) are complex. f On Poisson's hypothesis as to 
the relation between the elastic constants (X = p, h 2 = |-P), the roots of (62) are all 
real, viz., they are 

so that 

K jh-= -} v 7 (3 + v/3) = l e 087664 . . . ; 

this will usually be taken as the standard case for purposes of numerical illustration. 
In analogy with (28), it will be convenient to write 

where c denotes the wave-slowness of the Rayleigh waves. The corresponding 
wave-velocity is 

According as we suppose \ = co , or X = fx, this is *9553 times, or '9194 times, the 
velocity of propagation of plane transverse waves in an unlimited solid. 

The further properties of free Rayleigh waves are contained in the formulae (61) 

* 'Proc. Lond. Math. Soc./ vol. 17 (1885); 'Scientific Papers, 3 vol. 2, p. 441. 

t Of. Rayleigh (Joe. cit), where it is also shown (virtually) that they are roots of / (£), not of F (£), if 
a, /3 be chosen so as to have their real parts positive. 
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and (30). We merely note, for purposes of reference, that if in (33) we put ^ — i /c, 
and accordingly, from (61), 

A = (2^ - F) 0, B = ± 2iKa x G ...... (64), 

we obtain by superposition a system of standing waves in which 

u Q = — 2k {2k 2 — W — 2a 1 /3 l ) G sin kx . e^, t' = 2l<?a x G cos k# . e'?* . (65), 

The theory here recapitulated indicates the method to be pursued in treating the 
definite integrals of Art. 5. We fix our attention, in the first instance, on their 
"principal values," in Cauchy's sense, and afterwards superpose such a system of free 
Rayleigh waves as will make the final result consist solely of waves travelling 
outwards from the origin of disturbance. 

It may be remarked that an alternative procedure is possible, in which even 
temporary indeterminateness is avoided. This consists in inserting in the equations 
of motion (23) frictional terms proportional to the velocities, and finally making the 
coefficients of these terms vanish. This method has some advantages, especially as 
regards the positions of the "singular points" to be referred to. The chief problem 
of this paper was, in fact, first worked through in this manner ; but as the method 
seemed rather troublesome to expound as regards some points of detail, it was 
abandoned in favour of that explained above. 

7. The most important case, and the one here chiefly considered, is that ol a 
concentrated vertical force applied to the surface, to which the formulae (52) relate. 
The case of a horizontal force, expressed by the formulae (55), could be treated in an 
exactly similar manner. 

Since ti is evidently an odd, and v an even, function of 'x 9 it will be sufficient to 
take the case of x positive. 

As regards the horizontal * displacement u , we consider the integral 

J w ^ ~ J (2C - hj - v(£* - w) <y\p '- k*)£» "" * * [ h 

taken round a suitable contour in the plane of the complex variable £, = g + ^V 
If this contour does not include either "poles" (± k, 0), or "branch-points" 
(dt h, 0)> (± \ 0) of the function to be integrated, the result will be zero. 

A convenient contour for our purpose is a rectangle, one side of which consists of 
the axis of £ except for small semicircular indentations surrounding the singular 
points specified, whilst the remaining sides are at an infinite distance on the side 77 >0. 
It is easily seen that the parts of the integral due to these infinitely distant sides 
will vanish of themselves. If we adopt for the radicals \/(£' 2 — h 2 ) and \Z(C 2 — & 2 ), 

* The sense in which the terms " horizontal ;; and "vertical" are used is indicated in the second, 
sentence of the Introduction. 
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, « «a «m »™* /-^ 






Fig. 1. 



at points of the axis of £ the consistent system of values indicated in fig. 1 3 # we 
find, for the various parts of the first-mentioned side,! 



J — 00 e 



fr 



F(-k) 



" —a; J —k 






f *(£)dC=.iP 



p f (2f a - P^- 2a/3)e iix d ^ ■ k (2k 2 - g - g a ^) e 



sk# 



•A 



F(0 



F(/c) 



where the terms with F' ( — k) and F'(k:) in the denominator are due to the small 
semicircles about the points (± k, 0). Equating the sum of these expressions to 
zero, we find, since F ; (-~ k) = — F'(k), 

3 r fig^-^j-^^ = _ 2 ^ H cog ^ 



F(0 



«u1mm 



■-ft 



2 (* - £*_- 2*§ 2| 3 _- F + 2«/3 



F(f) 



/(O 



> &*** dg 



= - 2iM co* * - W {' £&£=$$^ .... (67), 

* The function under the integral sign in (66) is uniquely determined (by continuity) within and on 
the contour when once the values of the radicals J(£ 2 - h 2 ) and J(( 2 - ¥) at some one point are 
assigned. The convention implied in the text is that the radicals are both positive at the point ( + go , 0). 

It will be noticed that over the portion of the axis of £ between - h and - h the function in (66) 
differs from that involved in the value of uq as given by (52). This is allowed for in the second member 
of (67). Corrections, or rather adjustments, of this kind occur repeatedly in the transformations of this 
paper, 

t The symbol 13 is used to distinguish the "principal value" of an integral (with respect to a real 
variable) to which it is prefixed. 
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where 



H = 



k (2k 2 — k* — 2« 1 /3 1 ) 

F («) 



a • 



(68), 



a numerical quantity depending only on the ratio A : /a. 
To examine the value of v we take the integral 



*(£K = 



k^j^- h^e^clC 



(2£9 - #>)* _ 4 v/(r - 4' 2 ) ^(£ 3 - F) ^ 



. . (69) 



round the same contour. Integrating along the axis of £ we find 



'%«)««=» r-j«--«-» 



& 2 a t 



•GO 



■00 



F(f) 



F(-„) 



-e 



■J./C0C 






*(£)<*£ = 






.00 



-A 



*(£)<# = E 



/•00 



-A 



iPae** 36 d£ . A 2 a, 

5 — <?77« 1- 

F (£) w F (*) 



ITT ^.-t^: e lKX y 



and thence by addition, since the terms due to the infinitely distant parts of the 
contour vanish as before, 



where 



ff_ 



% ?ae^ x di; __ 



00 



F(0 



2iirK. cos kx + |P 



*f —CO 



+C{ 



1 1 



F(f) ' /(f) J 



k 2 ae^ x dg 



= — 2wK cos /c# + 2|| 






+ 2A; J r(af 7^ - (70), 



K 



F(k) 



(71). 



Hence if to the principal values of the expressions in (52) we add the system of 
free Rayleigh waves, 



u. 



o 



_ „• Q 



i — H sin /cr, 



-v, 



//, 







„• Q 



.ix 



K cos /ex . . ... (72), 
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which is evidently of the type (65), we obtain, on inserting the time-factor, 



u 










t' = 






.00 



% 



(2P 






/c 2 ) 



7T/x J ft (2P - A 2 ) 4 + 1 6£* (^ - h*) (& 3 - I 2 ) 



» b e 



(74). 



This is for x positive ; the corresponding results for x negative would be obtained by 
changing the sign of x in the exponentials, and reversing the sign of u . 

The solution thus found is made up of waves travelling outwards, right and left, 
from the origin, and so satisfies all the conditions of the question. 

The first term in u gives, on each side, a train of waves travelling unchanged with 
the velocity c~ l . The second term gives an aggregate of waves travelling with 
velocities ranging from 6" 1 to of 1 . As x is increased, this term diminishes indefinitely, 
owing to the more and more rapid fluctuations in the value of e^ iC . 

On the other hand, the part of v which corresponds to the first term of n Q remains 
embedded in the first definite integral in (74). To disentangle it we must have 
recourse to another treatment of the integral jM* (£) d£. One way of doing this is to take 
the integral round the pair of contours shown in fig. 2, where a consistent scheme of 



iZ 



^tt^^*m&&Mx*dtis*rrrza~ 



,f,tf. 



ivdv^r)*') 



i V(k + if) 



A 



£<t '-x-i&ifV. -m\2*^£ ^^^JKSJ^Si^W^Ca^SJtf^rtESifiitPM- J- 



"Sjjfr ^i. m - i, . hi j MM .-ffi fr 



>i/<h 2 i-y 2 ) 



iVilt^rf) 



% 




cc,fl 



a,-j3 -cc,-fi 



C>v ^ Lj 



-X. J- pL » Jmm a 



values to be attributed to the radicals \Z(t? — h 2 ) and \/(£ 2 ~~ &) is indicated. For 
the only parts of the left-hand contour which need be taken into account we find 
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"> (o <« = » r. %r - " ph.) e "' 
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F(0 
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Similarly, in the right-hand contour, 






ro 



i \/(^ 2 4- V 2 } k 2 e ^i drj 



. x (2f + F) 2 - 4t? 2 ^/(F + 7? 2 ) y(/c 2 + 7? 2 ) 



/•OO 







^P <T> dt = m [" — ?-—-£ — in ~*-~i e** 



We infer, by addition, 

i r £«£# = 27rK sin KX + 2 p ^z^f 



— 00 



F(6 



F(0 



J 7s 



1 



1 



F(f) /(f) 



lc 2 ae~ i ^d^ — 2 



r oo 



W sJ{W + 7? 2 ) C~^ (ZlJ 



o (2f + F) 2 - V v/(/^ + V 2 ) ^/(F + f ) 



2 T (75). 



If we multiply this by — Q/27r/x, and add in the term due to the free Rayleigh 
waves represented by (72), we obtain, as an equivalent form of (74), 



v 
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•h 



p /(#s - g*) e^'-M c? g 



TTfl 

4iQ f * 

7T/X Jfc 



o (2f 2 - F) 2 + 4£ 2 V /(F - £ 2 ) v/(F - £ 2 ) 



(2£ 2 - F)* +16^ (£ 2 - F) (F - £ 2 ) 



+ .Q. e ** 



»oo 



; 2 y/ (F + ?? 2 ) e-"* cfy 



tt/F Jo (2ij 8 + F) 2 - 4r, 2 v/(F + t? 2 ) y(F + r/) * 



(76). 



# 



It is evident that all terms after the first diminish indefinitely as x is increased. 

* From this we can deduce, by the same method as in Art. 4, an expression for the mean rate W at 
which a vertical pressure Q cospt does work in generating waves, viz., 



W = 






W ^(h 2 -J±d£ 



(2f - Wf + 4£ 2 J(h* - ? ) v/(^ 2 - £*) 






VOL. CCI1I. — A, 



W f* Ff ( £* - W) y/ (A 2 - £ 2 ) d£ 

D 
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If in (73) and (76) we regard only the terms which are sensible at a great distance 
from the origin, we have, for x positive, 

u =- Q -Re*^, v^^i-^ Ke^~-> . . ■ . . . (77) ; 

and similarly for x negative we should find 

u = Q He*<' +€ '> , v = - iQ K^ ( ' + "» . . . . . (78). 

These formulae represent a system of free Rayleigh waves, except for the 

discontinuity at the origin, where the extraneous force is applied. The vibrations 

are elliptic, with horizontal and vertical axes in the ratio of the two numbers 

H and K, which are defined by (68) and (71), respectively. To calculate these, we 

have, since F (k) = 0, 

/(*) = 2 (2k 2 - ¥f = Sa^K 2 , 
and therefore 

IT - ] l ( 2k * _ F ) 8 K - 2Fa i ^ - k ^ (7c>\ 

11 ^ _*F (*)/(*)' -F (*)/(*) • • ' • ' V /y * 

where, by differentiation of (62), 

-F(K)/( K ) = 16^{l-(6-4^)^ + 6(l-^)^} . . (80). 

In the case of incompressibility I find 

H = -05921, K = '10890; 

whilst on Poisson's hypothesis 

H = 12500, K = 18349, 

so that the amplitudes are, for the same value of jjl and for the same applied force, 
about double what they are in the case of incompressibility. 

A similar treatment applies to the formulae (55), which represent the effect of a 
concentrated horizontal force ~Pe ipt . Taking account only of the more important terms, 
I find, for x positive, 

a/ — _ ?i HVp**-^ 11 — — TT'p*P(t-cx) /Q1\ 



and, for a; negative, 

w = -^?HV* ( ' +ca!) , v =--KV^ +tt) . .... (82), 
/* . ' A l 

where 

tt - - **& - ^^ ( 2 ^nAT 

~~ F (k) ~ - ¥\k) /(«)" 



~,._ _ k (2k 9 -ft 8 - 2a, ft) _ A 8 (2 k 2 - A: 8 ) 3 



t I 00}« 
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The ratio of H 7 to K A is, of course, equal to that of H to K; K' is, moreover, 
identical with H, in conformity with the principle of reciprocity already referred to. 
It appears, therefore, from the numerical values of H, K above given, that for X = oo 



and for X = jx 



H 7 .= -03219, K' = -05921; 
B 7 = -08516, K' == -12500. 



Again, in the case of the internal source (56) I find, for large positive values of x, 



it 



o 



- 8/fHV^ ( 'r"), v & =. 8^KV ai ^^ a l .... (84), 



and, for large negative values, 



u 



o 



StcWe^e 1 *^^, v Q = 8iKK f e- a > f e ip(t+€x) ..... (85). 



The factor e~ aif indicates how the surface effect (at a sufficient distance) varies 
with the depth of the source. 

8. If in any of the preceding cases we wish to examine more closely the nature and 
magnitude of the residual disturbance, so far as it is manifested at the surface, it is 
more convenient to use the system of contours shown in fig. 3. With this system we 
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Fig* 3. 



can so adjust matters, that the radicals ^/(t? — K z ) and y/{t% — P) shall assume in all 
parts of the axis of £ exactly the values a, ]3 with, which we are concerned in formulae 
such as (52). It is convenient, for brevity, to denote by ± a, ft' the values assumed 
by the same . radicals on the two sides of the lines g.= .—:h 9 and by a", db/3" their 
values on the two sides of the line £ = — k, these values being- supposed determined 

D 2 
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In accordance with the requirements of continuity. Thus, with the allocation shown 
in the figure, we shall have, for small values of rj, 
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^/(2h v )e-^, /¥ = i v/(A* - ¥) 
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it 
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ft 
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approximately. 

Taking the integral 
the arrows, we find 




round the several contours, in the directions shown by 
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where, in the first integral, £ = — h + ^, and, in the second, £ = — h + ^. 
The integral (69), taken round the same contours, gives 
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(88), 



on the same understanding. 

The definite integrals in these results can all be expanded in asymptotic forms by 
means of the formula 



.CO 







,* x (v) e-*d n - 5L& x (o) + -W/? + ^ 4^ + • ' ' ' (89) ; 



and when fix, and therefore also kx 9 i& sufficiently large, the first terms in the 
expansions will give an adequate approximation, 
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Thus, taking account of (86), the last members of (87) and (88) are equivalent to 



2wH cos kx + 2 v/(2tt) a/ ( 1 



ky {kx)* 
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& s / (fcr) 1 



respectively. Substituting in (52), and adding in the system (72) as before, we have, 
for large positive values of x> 



u, = - Q He^^ + $ a/ 2 a/(i - ~ 



2\ gi (pt - i'X - \tt) 






*0 



'^(pt-hx-l-n) 



I " / \/ I 7 9/ / 7 V 



Q 



/ 2 m ^ ^'-fo-fr> 

T 2^V f ' (jfe» - 2A»)* * (for)* + *c. . . . (Jl j. 

The first terms in these expressions have already been interpreted. The residual 
disturbance constitutes a sort of fringe to the cylindrical elastic waves which are 
propagated into the interior of the solid, and consists of two parts. In one of these 
the wave- velocity p/k, or 6" 1 , is that of equivoluminal waves; the vibrations (at the 
surface) are elliptic, the ratio of the vertical to the horizontal diameter of the orbit 
being 2^/(1 — h?/k*) 9 or T633 for X = p. The remaining part has the wave-velocity 
p/h, or a" 1 , of irrotational waves ; the surface vibrations which it represents 
are rectilinear, the ratio of the vertical to the horizontal amplitude being 
(F — 2h^)/2h (k 2 — Wf y or *3535 for \=*[l. With increasing distance x the 
amplitude of each part diminishes as cc"" 1 , whereas in an unlimited solid the law 
is x~" i > as appears from (42)* 

Similar results will obviously hold in the case of the other problems considered in 

9. It has been assumed, up to this stage, that the primary disturbance varies as a 
simple-harmonic function of the time. It is proposed now to generalize the law 
of variation, and in particular to examine the effect of a single impulse of short 
duration. From this the general case can be inferred by superposition. 
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It is to be noticed, in all our formula, that if we write 

£ z= pO, h=£>ci, k = pb, k = pc, 

the symbol p which determines the frequency will disappear, except in the 
exponentials ; this greatly facilitates the desired generalization by means of 
Fourier's theorem. Thus, in the case of a concentrated vertical pressure Q (t) acting 
on the surface, the formulae (73) and (74) lead to 



v 



o 



Q («_,«)_ .^^^^Z^^^^^^^^ . Q(t-0x)d$.(92), 



H _ 



7r/Ja(2^--& 2 ) 4 + 16#*(# 2 -~a' 2 )(& 3 --^) ^ K ; 



1 «f & 9 v^ 2 - a 8 ) 



7T/x ^J* (2# 3 - & 2 ) 2 - 40* v /(<9 3 - a 2 ) v/(0 2 - 6 2 ) 



Q (i - &u) d9 . (93). 



The definite integrals represent aggregates of waves, of the same general type, 
travelling with slownesses ranging from a to 6, and from b to oo , respectively. 

If we suppose that Q (t) vanishes for all but small values of t, it appears from 
(92) that the horizontal disturbance at a distance x begins (as we should expect) 
after a time ax, which is the time a wave of expansion would take to travel the 
distance ; it lasts till a time bx, which is the time distortional waves would take to 
travel the distance ; and then, for a while, ceases. # Finally, about the timecse, comes 
a solitary wave of short duration (the same as that of the primary impulse) represented 
by the first term of (92). This wave is of unchanging type, whereas the duration 
of the preliminary disturbance varies directly as x,- and its amplitude (as will be seen 
immediately) varies inversely as x. 

If we put 

V^f ZZZ 1 \^ \0 J (a/0 • « . • » * • » » • \ /? 

the integration extending over the short range for which Q is sensible, the 
preliminary horizontal disturbance will be given by 



u 



_ ^H TT / ^ \ (Q^\ 

o - -LZhZ f U r™ ] ......... \vo), 

provided 

TT (f)\ - - h - 9 ( 2 ^^ ^h/i * ~ a ") ^i h l^l) (M) 

K ' ~~ . (20» - 6 2 )* + 1 6 &J& - a*) {W - 0*) ' ' ' ' V h 

where a <9 <b. The following table gives the values of U (0) for a series of values 
of 0/a, ori the hypothesis of X ==>, or b/a = 1*7321, 

* This temporary cessation of the horizontal motion is special to the case of a normal impulse. If 
the impulse be tangential, the contrast between the horizontal and vertical motions, in this respect, is 
reversed. 
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0/a. 


U(0). 


0/a. 


U(0). 


0/a. 


U(0). 


0/a. 


U(0). 


1-000 





1-025 


+ •62777 


1-10 


+ •22789 


1-550 


- -15122 


1-001 


+ -31247 


1-030 


+ •59351 


1-15 


+ ' 


10295 


1' 


600 


- -15842 


1-002 


+ -42080 


1-035 


+ -55806 


1-20 


+ ' 


02722 


1' 


625 


- -15927 


1-003 


+ -49148 


1-040 


+ -52308 


1-25 


— * 


02311 


1* 


650 


- -15681 


1-004 


+ •54191 


1*050 


+ -45741 


1-30 
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1 


•675 


- -14845 


1-005 


+ -57926 


1-060 


+ •39889 
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— ■ 
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- -12795 


1-010 


+ •66493 
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+ •34746 


1-40 


— < 


10771 
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•725 


- -07021 


1-015 


+ •67536 


1-080 


+ -30238 


1-45 
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•12527 


b/a 





1-020 


+ •65744 


1-090 


+ -26279 


1*50 


- -13975 




" 



The function has a maximum value + '67643 when 0/a = 1*01368 ; it changes sign 

when 0/a = 1*22474 ; and it has a minimum value — '159319 when 0/a = 1*62076. # 

The graph of this function is shown in the upper part of fig. 4. If the scales be 



a 









Fie- 4 

properly chosen, the curve will represent the variation of u with t 9 during the 
" preliminary" disturbance, at any assigned point x. For this purpose the horizontal 
scale must vary directly, and the vertical scale inversely, as x. 

* The calculations were made almost entirely by Mr. H. J. Woodall, to whom I am much indebted. 



24 



PROFESSOR HORACE LAMB ON THE PROPAGATION OF 



The interpretation of the expression (93) for the vertical displacement r is not 
quite so simple. For a given value of x, the most important part is that corre- 
sponding to t = ex, or = c, nearly, when the integrand in the second term changes 
sign by passing through infinity. This is the epoch of the main shock ; the minor 
disturbance which sets in when t = ax leads up continuously to this, and only dies 
out gradually after it. 

As a first step we may tabulate the function V (0) defined by 



V (0) = - 
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- -00193 


1-650 


- -44110 


1-004 


- -37630 


1*050 


- -25142 


1-30 


- -01508 


1-675 


- -52493 


1-005 


- -40039 


1-060 


- -20681 


1-35 


- -03796 , 


1-700 


- -63087 


1-010 


- -44907 


1-070 


- -16932 


1-40 


- -06941 


1-725 


- -76935 


X " Ul 


- "44543 


1-080 


- -13795 


1-45 


- -10989 


b/a 


- -81649 


1-020 


- -42324 


1-090 


-"* X 1 X i o 


1-50 


- "16137 


**"■"*""'"" 


«*^r—..st 



Ola. 



b/a 
1-75 
1-80 

X * oO 

6*M 



V(0). 



0-81649 
1-39031 
2-98197 
8-65843 

so 



Oja. 



1-90 
1-95 
2-00 
2*05 



V(0). 



+ 20-38685 
+ 5-42335 
+ 3*31759 

4- 2-46398 



01a, 



2*10 

J ' 10 

2-25 



V(0). 



+ 1-99591 
+ 1-69743 

+ 1 '48891 
+ 1*33404 



Ola. 



2-5 

3-0 

4-0 

10-0 



Y(0). 



+ -91464 
+ -60196 
+ -38179 
+ • 13292 



The function has a minimum value — '45120 when 0/a = 1*01170, and a zero 
maximum when 0/a = 1*22474; it changes from — oo to + oo when 0/h = 1*08767 ? 
or 0/a = 1*88389. # Its graph is shown in the lower part of fig. 4, and also (on a 
smaller scale, so as to bring in a greater range of 0) in fig. 5, 

It is postulated that the function Q (t) is sensible only for values of t lying within 
a short range on each side of ; the function Q (t — 0x) will therefore be sensible 
only for values of in the neighbourhood of t/x. We will suppose that for given 
values of x and t its graph (as a function of 0) has some such form as that of the 



* As In the case of U (0), the calculations are due chiefly to Mr, Woodall.' 
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dotted curve in fig. 5. If x be constant, the effect of increasing t will be to cause 
this graph to travel uniformly from left to right ; and if we imagine that in each of 



W) 



a 




c 








Fig. 5. 

its positions the integral of the product of the ordinates of the two curves is taken, 
we get a mental picture of the variation of v Q as a function of t, on a certain scale. 

For the greater part of the range of t, the integral will be approximately- 
proportional to the ordinates of the curve Y (0), viz., we shall have 



v. 







_y__ y (t 

TT/lbx * \X, 



(98), 



in analogy with (95), But for a short range of t, in the neighbourhood or cx 9 the 
statement must be modified, the dotted curve being then in the neighbourhood of the 
vertical asymptote of the function Y (0). Since the principal value of the integral is 
to be taken, it is evident that as t approaches the critical epoch and passes it, v will 
sink to a relatively low minimum, and then passing through zero will attain a 
correspondingly high maximum, after which it will decrease asymptotically to zero, 
the later stages coming again under the formula (98). 

Although the above argument gives perhaps the best view of the whole course of 
the disturbance, we are not dependent upon it for a knowledge of what takes place 
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about the critical epoch ex. We may proceed, instead, by generalizing the expres- 
sions (77). This introduces, in addition to the given function Q (£), whose Fourier 
expression is 



1 



,.00 



Q (t) = 1. 1 dp 

TT JO "" 



/»QO 



V —00 



Q (X) Qosp(t — A) <i\ 



ft 3 



• (99), 



the related function 



/»CO 



Q v (0 



/»00 



77 JO 



cZp Q (X) sin p (t — \)d\ . . . . . (100) * 



•^ ~oo 



viz., we have 



v 



o 



H 



Q (t — etc) + & c « 



v 







Q % (£ — ex) -f- &c. . . (101) 8 



It does not appear that the connection between the functions Q (t) and Q N (t) has 
been specially studied, although it presents itself in more than one department of 
mathematical physics. The following cases may be noted as of interest from our 
present point of view : 



\ql [tj — 



Q(*) = 



Q T 



77 t -f- T" 



2 ' 



Q(t) 



t 






# 9 * « * » 



(102); 



/ twwAieiew f*| 



Q v (0 



f »JL. i" 



iva 



, . . » ( 1 U o ) I 



Q(0 



_Q 



*a 



for £ a < t 



2 



2r 

for £ 2 > t 3 . 



f 



Q % (* 



log 



4?TT ° \t 



/ '£, — j— 7"\ 



T, 



(104), 



It is evident, generally, that if Q be an odd function, Q v will be an even function, 
and vice versa. 

The values of u and % as given by (101), are represented graphically in fig. 6, 
for the case where Q (t) and Q x (t) have the forms given in (102), # Moreover, writing 



we have 



HQ/2 tt/xt = /, KQ/27T/JIT = g, t — ex = t tan ^, 

u — — (1 + cos 2^) ./, ^o = s i n ^X-9 - 



. (105); 



the orbit of a surface-particle is therefore an ellipse with horizontal and vertical semi- 
axes f and g. And if from the equilibrium position O we project any other position 
P of the particle on to a vertical straight line, the law of P's motion is that the 
projection (R) describes this line with constant velocity. See fig. 7, where the 
positive direction of y is supposed to be downwards, 

* The relation between the scales of the ordinates in the graphs of u and v depends upon the ratio of 
the elastic constants A, /x. The figures are constructed on the hypothesis of \ — /a. . 
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Fig. 7. 



Fig. 6. 



A similar treatment would apply to the formulae (81), and (with some modification) 
to (84). 

It remains to justify these approximations by showing that the residual disturbance 
tends with increasing x to the limit 0. For this purpose we have recourse to the 
formulae of Art. 8. As a sufficient example, take the second term in the last 
member of (88). If we multiply by eM\ take the real part, and substitute r/ = p(j>> 
k = pb 9 the corresponding term in the value of % as given by (52), assumes the form # 

Q coap (t - bx) ( F (<£) e~*°* dcj> + Q. sin p (t - bx) f f(<J>) e~ p ** <#, 

[A J /X Jo 

where the functions F (<£) and /*(<£), which do not involve p, are of the order ^>~ 1 when 
<j> is large. If we generalize this expression by Fourier's Theorem (s^e equation (99)), 

we obtain, in the case of an impulse Q of short duration, 



-Q [°°F U) <te 

7T/X Jo 



TTfJL 



t*oo 



e 



—x$p 







cos p (t — 6x*) rfp + — ,/ (</>) ^ <?"* sin p(t — bx) dp 

7TUL J Jo 



7T/X Jo 



/•* 







FW 



cc</> r/<£ 






, 1106). 



* The symbols <f>, F, / are here used temporarily in new senses. 
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For any particular phase of the motion, t varies as x, and the expression (106) 
therefore varies inversely as x. This confirms, so far, our previous results (95) and 
(98). Hence with increasing distance from the origin the disturbance tends to the 
limiting form represented by (101). 

Before leaving this part of the subject, it is to be remarked that the peculiar 
protracted character of the minor tremor which we have found to precede and follow 
the main shock is to some extent special to the two-dimensional form of the question. 
It is connected with the fact, dwelt upon by the author in a recent paper, # that even 
in an unlimited medium a solitary cylindrical wave, whether of the irrotational or 
equivoluminal kind, is not sharply defined in the rear, as it is in front, but is prolonged 
in the form of a " tail." In the three-dimensional problems, to which we are about 
to proceed, this cause operates in another way. The internal waves are now 
spherical instead of cylindrical, and so far there is no reason to expect a protraction 
of a disturbance which in its origin was of finite duration. But at the surface they 
manifest themselves as annular waves, and accordingly we shall find clear indications 
of the peculiarity of two-dimensional propagation to which reference has been made. 
On the whole, however, it appears that the epochs of arrival of irrotational and 
equivoluminal waves are relatively more clearly marked and isolated than in the two- 
dimensional cases, 



PART II 



Three-Dimensional Problems. 



10. Assuming symmetry about the axis of z 9 we write 

= x /(x~ + y 2 ), u = — q, v = y~q .... (107). 



t»T 



so that q denotes displacement perpendicular to that axis. 

A typical solution of the elastic equations, convenient for our purposes, is derived 
at once from Art. 3, if we imagine an infinite number of two-dimensional vibration- 
types of the kind specified by (25) and (30) to be arranged uniformly in all azimuths 
about the axis of z, and take the mean. In this way we obtain from (33), with the 
necessary change of notation, 



9o = (ifjA - #B) . — ['e***™" cos cm da = — (fA + i£B) J\ (£W) 

77* J 



1 ' 

IT . 



>. (108), 



7T 



w {) = (— «A. — if B) . — . c^ C0S0J do) = — (a A + if B) J (fw) 

_ o 



Cited on p. 3/ post. 
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Also, from (40), for the corresponding stresses at the plane 2 = 0, we have 

O„] = M {2£«A + i (2f* - #>) B} J, (&r) ] 

> (109). 

[pj = /* { (2f» - F) A + 2#0B} J (£*) J 

Although the above derivation is sufficient for our purpose, it may be worth while 
to give the direct investigation, # starting from the equations 

3^?/ 3A 3 V 3 A 

P aF = (X + /0 9-- + /*V%, p ~- = (X + /*) -^. + ,,V*o, 

p^ ; = (X + ^)|| + /xV% . . (110), 

where 

A = ^ + -^ .4. 8 '~ (in) 



# 



In the case of simple-harmonic motion (e ipl ) these are satisfied by 

« = |*+V, v = d ^ + ,/, „ = &. + */. . . . (112), 

dx 3y 02; x /? 

provided 

(V* -f A 2 ) <£ = (J 13), 

and 



(V* + F) w' = 0, (V 2 + P) xf = 0, (V 2 + A 2 ) u/ = 



-\ 



3tt' , 3o' , 3 jo' „ f ■ • • ( 114 ): 

-4— — r— __ (J 

dx dy dz 
where /i 2 , W are defined as before by (28). A particular solution of (114) is 

3^3;; 3t/3z 3^ v /x v 7 

provided 

(V 2 + F) x = . , (116). 

On the hypothesis of symmetry about Oz we have 



|3 



3 2 1 3 3 

czar 7aT dm oz z ; 



and the formulae (112), (115) are equivalent to 



? 



? = |i + 3!x, w 3* + ^ + ^ (U8) , 

3ct 3^3;; 82; 



* (?A ' Proc, Lond. Math. Soc./ vol. 34, p. 276, for the corresponding statical investigation, 
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If we take,, as the typical solution of (113) and (116), 

<f> = Ar w J (£w), x = B«r* J (&r) ..... (119), 

where a, /3 have the same meanings and are subject to the same convention as in 
Art. 3, we have, from (118), 

q = (- fAe-« + fiSBeT*) J x (&r) 
w = ( — aAc"""* + £ 3 B<r p *) J (£rar) f 

and thence for the stresses in the plane z = 



a . . » . ( L Zl U I , 



-PevrJO — 


fl 


dq JL div 
Jbz ' 3rar_ 





-P*sj0 = 


X 


1 + 2/x -~ 

dz 


r ; 







= /i {2faA-(2^-F)fB}J 1 (^) 
= /x {(2£ a - F) A - 2f s /8B} J (&r) 



^ . (121). 



y 



The formulae differ from (108) and (109) only in the substitution of i£B for B. The 
notation of (119) is adopted as the basis of the subsequent calculations. 
If we are to assume, in place of (1 1 9), 



<£ = A'e"J (£*), x = B'e*J (&r) • 



. I L £j Z: J . 



the corresponding forms of (120) and (121) would be obtained by affixing accents 
to A and B, and changing the signs of a and /3 where they occur explicitly. 

11. As in Art. 4, we begin by applying the preceding formulae to the solution of a 
known problem, viz., where a given periodic force acts at a point in an unlimited solid. 

Let us suppose, in the first place, that an extraneous force of amount Z . J (^W) e^ f \ 
per unit area, acts parallel to z on an infinitely thin stratum coincident with the 
plane z = 0. The formulae (119) will then apply for z > 0, and (122) for z < 0. The 
normal stress will be discontinuous, viz. : 



LiV>-+0 iP^jz 



Zj e iJ A ( £13X ) 



whilst p^ is continuous. Hence 



(2£ a - ¥■) (A - A') - 2P/3(B + B') 



Also, the continuity of q and ^ requires 



A - A - /3 (B + B') = 
a (A + A') - £ 2 (B - B') 







We infer 



.0. 



A = 



JLk 



Z 



"N 



> 



2a (A + A') — (2p - fc 4 ) (B — B') = 



^ 



* * 



&jw 



B- B 



>/ 



'P 



&V$ 



(123) : 



(124). 



s e 



(125). 



a e 9 



(126), 
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and therefore, for z > 0. 



4> 



2/c 



Z Z e~^ 

e ~ as J o (&*) > X = o ji- -V J o (&0 



y* 



ZifC LL kj 



. (127). 



To pass to the case of a concentrated force Re^, acting parallel to z at the origin, 
we have recourse to the formula (20), where we suppose jf(X) to vanish for all but 
infinitesimal values of X, and to become infinite for these in such a way that 



/»QO 



f (\) 2it\ dk = II. 

We therefore write Z = T&£ d£/27r, and integrate with respect to £ from to oo , # 
■We thus find, for z > 0, 






E 



i.wp^p 



/»00 

J 



it f tf f 



ATrp^pU (3 



J (&r)fd£. . (128).. 



which are equivalent, by (18), to 



* 



13 



3 e -ihr 



Arrp^p ' dz r 



R 



e 



■^ 



X : — 



4:irp z p 



/y* 



* \ JL <Li if J t 



This will be found to agree with the known solution of the problem, f If we retain 
only the terms which are most important at a great distance r, we find, from (118), 






R 



-, j 



w 



4:TT 



R 



Zm 



• ikr 



Z7aT 



1 



fX T 



g — ikr I 1 



Z % _ 



4tt I X -f" 2/x 7 



/]»d 



e 



ihr i _„ ^ _ g — Mr 



/i 4 



-3 



!► 



. (130). 



-/ 



Inserting the time- fact or, the radial displacement is 



R 



z 



zw + zsq _ 

r, Arr (X + 2/x) * r 2 



gip a — «?•) 



(131), 



and the transverse displacement in the meridian plane is 



7j$W 



zq __ R 

4z7Tfl 



3?L fJ ip(l ~br) 



» « • . ♦ (.LoziK 



Returning to the exact formulas (128), the expression for the velocity parallel to 
at the plane z = is found to be 



z 



cw {Re^' jt 



dt 



4rrpp 



/»Xi 



JQ 



2 N 



a + i) Jo (£*)£<*£ ....... (133), 



* A more rigorous procedure would be to suppose in the first instance that the force R is uniformly 
distributed over a circular area of radius a, using the formula (22). If in the end we make a "~ 5 we 
obtain the results in the text. 

t Stokes, <Caml>. Trans./ vol. 9 (1849); < Mathematical and Physical Papers/ vol 2, p. 278. 
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or, taking the real part, 



dw 
dt 



r rk 

L J 



f 



rh 



Airpp I ' o i/(k* 



>gy J (f CT ) d£ + | £ \/( /i3 ™ f*) Jo (6») rff I cos pt 

"" I ufc/X JLlIo XXI oXXX L/v » a » • • * • • • 



(134). 



The terms in cos pt remain finite when we put or = ; # and the mean rate W at 
which a force R cos p£ does work in generating waves is thus found to be 



W 



R* 



£ 3 df 



8rrpp 

R2 

24irp/> 



/» 



y/(k? — £*) JO V V ' 

. (2P + /i 3 ) = ^ (a 3 + 2& 3 ) . . 



* * 9 I JL %J %J I * 



a and b denoting as before the two elastic wave-slownesses. The result (135) can be 
deduced, as a particular case, from formulae given by Lord Kelvin, f 

12. Proceeding to the case of a semi-infinite solid occupying (say) the region 2> 0, 
we begin with the special distribution of surface-stress : 

[p**\o = Z . J (£sr), [P*ar] = ° • ..... .'(136). 



The coefficients A, B in (119) are now determined by 



(2f* - F) A - 2f /3B 



f 



whence 



2aA-(2f s -Jfc*):B = () 



a * * # 9 



• • (137), 



J 



xL — — 



«2f 



2 



W 



F(a 



/* 



15 



a 



F( 



# » 



e s « 



z 1 



(138), 



the function F (£) having the same meaning as in Art. 5. The corresponding 
surface-displacements are 



?o = 



u> — 



£ (2f 2 - F - 2«/3) T 
F(£) 

¥ y } • Jo (*») • - 



(f») 



"> 



ft 



^ 



9 s a e © 



139). 



-^ 



This result might have been deduced immediately from (51) in the manner indicated 
at the beginning of Art. 10. 

* The terms in sin pt become infinite. If the force R be distributed over a circular area, the awkwardness 
is avoided. A factor 

Ji m \ 8 

is thus introduced under the integral signs in the first line of (135), where a denotes (for the moment) the 
radius of the circle. Finally, we can make a infinitely small 
t 'Phil. Mag./ Aug. 1899, pp. 234, 235. 
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Ti we put Z = in (137) we get a system of free annular surface-waves, in which 

q = — k (2/c 3 — P — 2^/3^ . Jj (kzj) . Ce^t 



w = /c 2 a 1 . J (kzs) . CeM 



y . . . . (i4o). 



where k is the positive root of F(£) = 0, and a ]5 /3 X are the corresponding values of 
a, /3. These are of the nature of " standing " waves. 

To pass to the case of a concentrated vertical pressure He¥ at 0, # we put in 
accordance with (20), Z = — llgd£/2TT. and integrate from to oo.f The formulae 
(139) become 



90 ""2^ Jo F(f)~ ~ * J i^)^ | 



^o 



27r M JoF(f) J oV^)^ 



J 



. (141). 



Again, the case of an internal source of the type 

— ihr 



P, 



<j> — - -• ■ , x ^ ^ • * • * * "• - - • (142), 

where r denotes distance from the point (0, 0, /), can be solved by a process similar 
to that of Art. 5. First, superposing an equal source at (0, 0, —-/), distance from 
which is denoted by r\ we have 

P — ihr p — ihr' 

0= + -, -, X : =° (143); 

and therefore, by (18), in the neighbourhood of the plane z = 0, 

<f> = J (fa) &£ + - J (fa) £tf£ 

J Q OL J o OL 



This makes 



2 | ^c^Jfl^)^ ......... (144). 



■«> -<*/■ 



g = — 2 " Ji(fa) ij' z dij, w = . . . . . (145). 

■ ' J o « 



* This may be regarded as the kinetic analogue of Boussinesq's well-known statical problem. 

t It might appear at first sight that a simpler procedure would be possible, and that the effect of a 
pressure concentrated at a point might be inferred by superposing lines of pressure (through O) uniformly 
in all azimuths, and using the results of § 7. It is easily seen, however, that such a distribution of lines 
of pressure is equivalent to a pressure-intensity varying inversely as the distance (w) from O. This is 
not an adequate representation of a localized pressure, since it makes the total pressure on a circular 
area having its centre at O increase indefinitely with the radius of the circle, 
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anc 



1 



" P-- Jo = 0. O Jo = 2 /* 



-X) 







(2f 



r 



a 



e~«h {&)£«£ 



(146). 



The additions to (1.43) which are required in order to annul the stresses on the 
plane z ; = are accordingly found by writing 

- jfc* 



2/x • 






a 



c-* i-dg 



in (139), and then integrating with respect to f from to oo . In this way we 
obtain, finally, 



/•OO 



<7o 



4l o |^^'J 1 (^)r^ 



"1 



.^oo 



w 







2 



^ (2f 







(£) 



ki hr«,J a (fc)d£ 



y 



* * 



(147). 



In a similar manner, with the help of Art, II, we might calculate the effect of a 
periodic vertical force, acting at an internal point, 

13. For the sake of comparison with our previous two-dimensional formulae, it is 
convenient to write, from (2) and (6), 



Jo (H 



Ji(£*) 



% 



,'jj 



it J a 
1 



(e« 



w cosh u 



er^*°* h "")du 



TV - 



"CO 



Y 



(148). 







The formulae (141) are thus equivalent to 



<?o 



/♦OO 



2u> 



>oo 



1 cQgh w r / w 







■ :o 






_^P) pit- ix cosh h .7£ I 



w 







2tt 2 /x 



,-»oo 



/■•CO 



of-w 



«/ 







J . — oo 



7 2 £ 



^ 



(149). 



^ 



These results are closely comparable with (52), and our previous methods of treat- 
ment will apply. It is, however, unnecessary to go through, all the details of the 
work, since the definite integrals with respect to £ which appear in (149) can be 
derived from those in (52) by performing the operation — id/dx upon the latter, and 
then replacing x by tsf cosh u. 

Thus, from (67) and (70) we derive 



»■» J-C 



L ( 2 f . A . 2a P) e i^coBhv c ig _ 2tt/cH sin (kot cosh w) 



•Xi 



F(f) 
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f I ■ I'fe e** cosh u d£ = 2tt/cK sin ( K m cosh w) - 2F$ 
-*> " (&) 



£ 



=> p — it-Zu COSlli' 

F(0 



c# 



op f " fe (2? — « )"a -,| W cosh« ,/dr / 1 c i \ 



where H and K are the numerical quantities defined by (68) and (71). Substituting 
in (149) we have 



$<io = 



. H . Kj (k-ot) + -i? //-\ >/>\ "'■" l*i (& w ) tt 4" • (1^2). 



2/i, 



7r^t 



/« 



F (£)/(£) 



«K 



■r»0O 



:]f}w = ----'■- . K . J (fccr) 

A i / / 



jJL 



t* »i; f|) d » <« * 



"SfL"",?! 3 "l)o(^)^. • (153) : 



_ iWR, 

2^ J* F (£)/(£) 

where the notation of the various Bess el's Functions is as in Art. 2. 
Superposing the system of free waves in which 



% = 2 . M. J x (/cot). 



w 







2/T 



. K . J (kts) . (154), 



we obtain, finally, on inserting the time-factor 



2o 



^.H.D i( k^+^ . (155), 






2/X 7TjLt J£ 



,%)/(£) " D » (*»)«*#. . (156). 



Since these expressions are made up entirely of diverging waves, they constitute 
the complete solution of the problem where a periodic normal force Ke^ is applied to 
the surface at the origin. 

An alternative form of (156), which puts in evidence that part of the vertical 
disturbance which is most important at a great distance from the origin, is obtained 
from (75). Attending only to the " singular" term, we find 



'X 



W& -n 



i I ^ p >" e^ooshH d g _ _ 2i7r/cK _ cos ^ cogh u j + &c _ 



• * 



(157), 



and therefore, from (149). 
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^Wq = - ■■■ . K . K (kzj) + &c. ....... (158). 



Adding in the system (154) we have altogether 
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Hence, by (7), we have, at a great distance w, 

q =- *^H. a/-? ..-. e -<P<-"'- W, | ( ; ^ kE K. a/- 2 -.^^-"- *>. (160). 
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This may be compared with (77). The vibrations are elliptic, with the same ratio of 
horizontal and vertical diameters as in the case of two dimensions ; but the ampli- 
tude diminishes with increasing distance according to the usual law zs~^ of annular 
divergence. 

In the same manner we obtain, in the case of an internal source of the type (142), 



q Q = — -~ rr f ~- 1 e~ a ' f D, (kot) e' 1 ' 1 -f- &c., 

iv () = ~- -A, / x a D n (/co5") e l + cs;c. 

u F'(k) o\ / « 
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where the factor e ~ a]J shows the effect of the depth of the source. 

The expressions for the residual disturbance might be derived from the formulae 
of Art. 8 by the same artifice. Without attempting to give the complete results, 
which would be somewhat complicated, it may be sufficient to ascertain their general 
form, and order of magnitude, when hrs and km are large. To take, for example, the 
parts due to the distortional waves, if we perform the operation — id/dx on the 
second terms of the unnumbered expressions which occur between equations (89) 
and (90), above, and then replace x by m cosh u, the more important part of the 
result in each case is 

e- ik ^ w " hu /(km cosh uf ll \ 

multiplied by a constant factor. This result is to be substituted for the definite 
integrals with respect to f which occur in (149); the corresponding terms in 
f/ and iv are therefore of the types 
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and 



(fe)Uo (cosh uf ' (kzsf 1 Jo (cosh uf 

respectively. By the method by which the asymptotic expansion (7) of the 
function D (£) is obtained, it may be shown, again, that these terms are ultimately 
comparable with 

e il)(t ~* bm) /(km)\ 

where the time-factor has been restored. In the same way, the terms in q Q and iv 
which correspond to the expansional w r aves are ultimately comparable with 

The attenuation with increasing distance is much more rapid than in the case of the 
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annular llayleigh waves, so that the latter ultimately predominate. # - It is also 
much more rapid than in the case of elastic waves diverging from a centre in an 
unlimited medium, where the amplitude varies inversely as the distance. 

14. The generalization of the preceding results, so as to apply to an arbitrary 
time-variation of the source, follows much the same course as in Art. 9. The full 
interpretation is however more difficult, so far at least as regards the minor tremors. 

The main part of the disturbance, in the case of a local vertical pressure applied to 
the surface, is obtained by generalizing the formulae (159). These may be written 

q = B ? A f e l ^- m ^ h ^du + &c, w = - lK lic | f e l i ){t - mco * hu) du + &c. (162). 

IT /X OZaT Jo ' TT [JL ot J 

Hence, corresponding to an arbitrary pressure R (£), we have 



h a 

% = — a- 

7T/X OCT 



ftt Y7~ <"\ »Q0 



R(£ — C7ST cosh u) du-\-&,c> w = — ----- j Ii v (£~ ctscosh u) du-^-ko,. (163), 

7T/X dt Jo 



where, in analogy with (100), 
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B, N (£) = — dp E (\) snip (t — \) d\ . . . . . (164). 
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The character of the function of t represented by the first definite integral in (163) 
has been examined by the author f for various simple forms of H (t), and a similar 
treatment applies to the second integral. For example, if we take 

R(*) = -^ 1 i> K X (0 = - / a4-T ...... (165), 

it is found, on putting 

t — . cot == r tan x,§ 

that for values of w large compared with r/c, and for moderate values of x> 
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11 (t — cot cosh w) du = ^ \/('~ I cos (i 77 ~~ ix) v / ( cos x) • (166)t» 



I ET (£ — cot cosh n) cZ/.t = — ■" : - \f ( ,T ) sin ({tt .■— ix) \/( cos x) • ( 16 ^)' 

Jo 2r v \Cot/ 

approximately. Substituting in (163), we have, ignoring the residual terms, 

20=— / sin (i 77 ~ tx) cosl X ~ 
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o 



# cos (^r - f x) cos f X 



>• . . ../ . . . (168). 



* Of. the footnote on p. 2 ante. 

t "On Wave-Propagation in Two Dimensions," ' Proc. Lond. Math. Soc.,' vol. 35, p. 141 (1902). 
I Of. Equation (36) of the paper cited. It may be noticed that the functions on the right hand of (166) 
and (167) are interchanged, with a change of sign, when we reverse the sign of x- 
§ The symbol x is no longer required in the sense of equations (115), &c. 
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where 
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The following numerical table is derived from one given on p. 155 of the paper 
referred to : — 
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The graphs of q () and iv as functions of t P in the neighbourhood of the critical 

epoch czar, are shown in fig. 8, which may be compared with fig. 6.f The corresponding 

orbit of a surface particle is traced in fig. 9, where the positive direction of % is 

downwards ; it may be derived by a homogeneous strain from a portion of the curve 

whose polar equation is 

r* = a 1 cos § (0 — • \ir\ 

The amplitude of this part of the disturbance diminishes, with increasing distance 
from the source, according to the law m~~K 

Complete expressions for the disturbance are obtained by generalizing (155) and 
(156). They may be written 
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where U (0) and V (6) are the functions defined and tabulated in Art. 9 



t See the footnote on p. 26 ante. 
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Fig. 9. 



Fig. 8. 

The method applied in that Article to obtain a general view of the whole progress 
of the vertical displacement at any point might be employed again here, the upper and 
lower curves in fig. 4 being combined with auxiliary movable graphs of 
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considered as functions of 9. In the case of a primary impulse of the type (165), both 
graphs would have somewhat the form of the lower curve in fig. 8, the functions being 
practically (except for a constant factor) of the type 



9 



- sin (|tt - |x) cos'x , 
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in the more important part of the range. Both graphs, if drawn to the scale of 
fig. 4 or 5, would be excessively contracted horizontally when we are concerned with 
values of w large compared with r/c . Owing to the compensation between positive 
and negative ordinates in the auxiliary graphs, it is plain that the disturbance 
expressed by the ^-integrals in (169) and (170) will be relatively very small except 
when t/m has values 9 for which the gradient of U (9) or V (6) is considerable. As 
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regards the horizontal displacement q () , the minor tremor will consist of a single 
to-and-fro oscillation about the epoch am, followed after an Interval by a somewhat 
similar oscillation about the epoch bm, with almost complete quiescence between. As 
regards the vertical displacement, there will be a to-and-fro oscillation, about the epoch 
am, then a period of comparative quiescence, and finally a gradually increasing 
negative displacement (with a slight irregularity at the epoch bm) leading up to the 
main shock, after which there Is a gradually decreasing positive displacement. 

The expression for the horizontal displacement q (] may be treated in a different 
manner. Transforming (169) we have 
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A rough sketch of the graph of tfU 7 ^) + U($) is easily made, and the function 
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is, in such a case as (165), one-signed, but its integral with respect to does not 
converge when the lower limit Is large and negative, The method therefore fails to 
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Fig, 10. 
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give us a convenient view of the progress of q as a function of t. The difficulty is 
due to the peculiarities of annular propagation to which reference has already been 
made. 

In fig. 10 an attempt, based on the former method, is made to represent (very 
roughly) the whole progress of the horizontal and vertical displacements due to a 
single impulse of the type (165) at a distance large compared with r/c. 

Summary. 

We may now briefly review the principal results of the foregoing investigation, so 
far as they may be expected to throw light on the propagation of seismic tremors 
over the surface of the earth. 

It has been necessary to idealize this problem in various ways in order to render it 
amenable to calculation. In the first place, the material is taken to be compact and 
homogeneous, to have, in fact, the properties of the " isotropic elastic solid " of 
theory. Moreover, the curvature of the surface is neglected. Again, instead of a 
disturbance originating at an internal point, we study chiefly the case of an impulse 
applied vertically to the surface. Under these conditions the disturbance spreads 
over the surface in the form of a symmetrical annular wave-system. The initial form 
of this system will depend on the history of the primitive impulse, but if this be of 
limited duration, the system gradually develops a characteristic form, marked by 
three salient features travelling with the velocities proper to irrotational, equi- 
voluminal, and Rayleigh waves, respectively. As the wave-system, thus established, 
passes any point of the surface, the horizontal displacement shows first of all a single 
well-marked oscillation followed by a period of comparative quiescence, and then 
another oscillation corresponding to the epoch of arrival of equivoluminal waves. 
The whole of this stage constitutes what we have called the " minor tremor" ; it is, 
of course, more and more protracted the greater the distance from the source, and 
its amplitude continually diminishes, not only absolutely but also relatively to that 
of the " main shock," which we identify with the arrival of the Rayleigh wave. It may 
be remarked that the history of the minor tremor depends chiefly on the time- 
integral of the primitive impulse ; the main shock, on the other hand, follows the 
time-scale of the primitive impulse, and is affected by every feature of the latter. >x< 

Similar statements apply to the vertical displacement, except that the minor 
tremor leads up more gradually to the main shock, being interrupted, however, by a 
sort of jerk at the epoch of arrival of equivoluminal waves. 

The history of the horizontal and vertical displacements, about the epoch of the 
main shock, in the case of a typical impulse of the type (165), is shown in fig. 8 ;. 

* Observational evidence in favour of the existence of the three critical epochs in an earthquake 
disturbance has been collected and discussed by E. D. Oldham, "On the Propagation of Earthquake 
Motion to Great Distances/ 5 'Phil. Trans., 5 A, 1900, vol. 194, p. 135. 
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whilst fig. 9 shows the corresponding orbit of a surface-particle. In fig, 10 a sketch 
is attempted of the whole progress of the disturbance. 

These results are of a fairly definite character, but they are based, as has been said. 
on purely ideal assumptions, and it remains to inquire how far they are likely to he- 
modified by the actual conditions of the earth. The substitution of an interned source 
for a surface impulse will clearly not affect the general character of the results at a 
distance great compared with the depth of the source, although differences of detail 
in the wave-profile at the critical epochs will occur, and we can no longer assume that 
the disturbance is the same in all vertical planes through the source. Again, the 
chief qualitative difference introduced by the curvature of the earth will be that 
the minor tremor, whose main features are evidently associated with the outcrop of 
spherical elastic waves at the surface, will be propagated directly through the earth, 
so that the first two epochs will (at distances comparable with the radius) be 
accelerated relatively to the main shock, # which, being due to the Hayleigh waves 
will travel, with the velocity proper to these, over the surface.]' 

It is a more difficult matter to estimate the nature and extent of the modifications 
produced by heterogeneity. It is, perhaps, possible to exaggerate these, for the 
qualitative effect of a gradual charge of elastic properties would not be serious, and 
even considerable discontinuities would have little influence if their scale were small 
compared with the wave-length J of the primitive impulse. A covering of loose 
material over the solid rock probably causes only local, though highly irregular, 
modifications, with some dissipation of energy. 

It must be acknowledged that our theoretical curves differ widely in two respects 
from the records of seismographs. In the first place, they show nothing corresponding 
to the long successions of to-and-fro vibrations which are characteristic of the latter. 
It would appear that such indications, so far as they are real and not instrumental, 
are to be ascribed to a succession of primitive shocks, in itself probable enough. 
Again, the theory gives vertical and horizontal movements of the same order of 
magnitude, and in the case of the Rayleigh waves, at all events, where a definite 
comparison can be made, the vertical amplitude is distinctly the greater. The 
observations, on the other hand, make out the vertical motion to be relatively small. 
The difficulty must occur on almost any conceivable theory, and appears indeed to be 
clearly recognised by seismologists, who are accordingly themselves disposed to 
question the competence of their instruments in this respect. 

* Of. R. D. Oldham, he. cit. 

t The theory of free Rayleigh waves on a spherical surface is known ; see Professor Bromwich, loc. cit. 

J This term is used in the same general sense in which in hydrodynamics we speak of the "length" of 
a solitary wave travelling along a canal. There is no question, in the present connection, of anything 
analogous to "oscillatory waves." 



